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G> ! Abstract 

S ■ Dynamics of a tracer particle subject to a constant driving force in a one- 

pH ■ dimensional lattice gas of hard-core particles whose transition rates are sym- 

o ■ 



O ■ metric. We show that the mean displacement of the driven tracer, XT{E,t) 



vj grows in time t as XT{E,t) = y/od, rather than the linear time dependence 

■ found for non-interacting (ghost) bath particles. The prefactor a is deter- 

mined implicitly, as the solution of a transcendental equation, for an arbitrary 
magnitude of the driving force and an arbitrary concentration of the lattice 
gas particles. In limiting cases the prefactor is obtained explicitly. Analytical 
predictions are seen to be in a good agreement with the results of numerical 
simulations. 
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I. INTRODUCTION 



Dynamic and equilibrium properties of lattice gases, i.e. systems involving randomly 
moving particles with hard-core interactions, have received much interest within the last 
several decades. A number of important theoretical results have been obtained for such 
systems revealing non-trivial, many-body behavior p|-p!2|. Lattice gas models often serve as 



microscopic models of complex physical phenomena. To name a few we mention dynamics of 
motor proteins |T3|,|T^, growth of interfaces p!5|-p!7|, traffic jams and queuing problems [0. A 



lattice-gas approach has been used for the derivation of Euler-type hydrodynamic equations, 
e.g. the Burgers equation [Il9ip0|. Another important example concerns the spreading of 



molecularly thin wetting films, where experimental studies pl|-|23| have evidenced surprising 
universal laws which recently have been explained in terms of a lattice gas model ||2^-|26 



We believe also that a robust microscopic, molecular approach to such physical phenomena 
as shear-induced ordering in colloidal suspensions |^ or stick-slip motion of mica planes 
separated by an ultrathin liquid layer [p^-|3T| could also begin with a description based on 



a lattice-gas picture. 

At present time two models are well-studied in the literature. In the first, the so-called 
asymmetric exclusion process, all particles in the system perform stochastic motion, con- 
strained by hard-core interactions, in the field of a constant driving force |l|-[9|,p^-|3^ . Here 



the velocity, diffusion constant and equilibrium configurations have been calculated exactly 
for different types of boundary conditions ( e.g. and references therein). In the 

second, no external force is present and all particles have symmetric transition rates. Re- 
markably, in such a situation, the motion of a labeled, tracer particle is non-diffusive in 



low dimensions. For example, the mean-square displacement X^{E = 0,t) of a tracer par- 
ticle (identical except for its observability to all other particles) in a one-dimensional (ID) 
symmetric lattice gas shows a sublinear growth with time p5|j3( 

Co /"^t^i/ 



xuE = o,t) = — ^(-)^/^ (1) 

1 - Co vr 

where Cq, < Cq < 1, denotes the mean (constant) concentration of vacant sites and 



the argument E = signifies that the external force is absent and that all particles have 
symmetric transition rates. Hence, in ID trajectories Xt{E = 0,t) of such a tracer particle 
are more compact than these of particles without the hard-core constraints. In 2D the 



mean-square displacement X^{E = 0,t) shows a linear dependence on time with additional 
logarithmic terms |Tl| , |37| ; in 3D it grows linearly in time with the diffusion constant being 
a nontrivial function of the particle concentration. 

In the present paper we focus on the less studied and less understood situation in which 
only one particle (the tracer) experiences the action of an external (constant) driving force, 
E and thus has asymmetric transition rates, while all other particles (bath particles) are 
not subject to this force and have symmetric transition rates. The tracer behavior in such 
a system with a vanishingly small driving force has been first examined in in which 
the question of the validity of the Einstein relation for the hard-core lattice gases has been 
addressed^. Such a model has been used [|38| in a numerical study of the gravity driven 



motion of a finite rigid rod in a "sea" of hard-core monomers. Another physical example 
corresponds to the situation in which a charged particle diffuses in a lattice gas of electrically 
neutral particles in the presence of a constant electric field E. The extreme case of infinitely 
strong electric fields {E = oo), which means that the tracer particle may move only in one 
direction, has been studied in [^. It has been shown, for example, that in ID systems the 



mean displacement, Xt{E = oo, t), of a charged tracer particle grows sub linearly with time, 

XTiE = oo,t) cx («oot)'/', (2) 



where Ooo is a constant |^. Eq.(2) indicates that hard-core interactions give rise to an 
effective friction. In the one-dimensional case this force is much stronger than the viscous 
friction for one particle and grows in proportion to the mean displacement of the tracer , i.e. 
as \/t, as this is a measure of the size of the compressed region preceding the tracer which 
hinders the ballistic motion of the driven particle. 



^This point will be discussed in more detail in Section IV. D. 
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Here we study motion of a driven tracer particle in a symmetric lattice gas in the general 
case of fields of arbitrary strength and at arbitrary concentrations of the lattice-gas particles. 
Focusing on one- dimensional situations in which the hindering effect of the lattice-gas (bath) 
particles on the tracer motion is most pronounced, we devise a mean-field-type theory which 
allows simple calculation of the mean displacement of the tracer particle as a function of 
time and other pertinent parameters. We find that XT{E,t) has the following dependence 

xtOM = (« ty/', (3) 

where the parameter a is a time-independent constant, which is a complicated function of 
field strength E, which determines the transition probabilities p and q, and concentration of 
the bath particles, Cp,. This constant is determined here for arbitrary values of E and Cp. 
Our analytical findings are in excellent agreement with the results of numerical simulations. 

The paper is structured as follows: In Section 2 we describe the model. In Section 3 
we present definitions and write down basic equations describing motion of particles. In 
Section 4 we determine explicitly the growth law for the mean displacement of the tracer 
particle, Eq.(3), and evaluate a closed transcendental equation for the parameter a. In sev- 
eral limiting cases the dependence of a on the pertinent parameters is explicitly obtained. 
Section 5 presents results of numerical simulations and comparison of these with our an- 
alytical predictions. Finally, In Section 6 we conclude with summary of our results and 
discussion. 

II. THE MODEL. 

The model is defined in the following way. Consider a one-dimensional regular lattice of 
unit spacing, infinite in both directions, the sites {X} of which are either singly occupied by 
identical particles or vacant. The particles are initially placed at random (constrained by the 
condition that double occupancy of sites is forbidden) with mean concentration Cp = 1 — Co, 
where Cp is the mean site occupancy, Cq being the mean site vacancy. The tracer particle 
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is put initially at the origin, i.e. at X — 0. A configuration of the system is characterized 
by an infinite set of (time-dependent) occupation variables {tx}, where tx = 1 if site X is 
occupied and rx = if site X is vacant. Consequently, the variable r]x = ^ — tx describes 
the probability that site X is vacant. 

The dynamics of the bath particles is symmetric: each particle waits a (random) expo- 
nentially distributed time with mean 1 and then attempts to jump, with equal probability 
(1/2) to the right or left neighboring site. The jump actually occurs if the chosen site is 
empty. The tracer particle motion is asymmetric: the tracer waits a random exponentially 
distributed time with mean 1 and then randomly selects a jump direction. It chooses right- 
hand (left-hand) adjacent site with probability p {q — 1 —p). The jump occurs if the selected 
site is vacant. If the asymmetry in tracer jump probabilities is due to an external electric 
field, E, one has the relation p/q = exp{(3E), where (3 is the inverse temperature. For sim- 
plicity we have set the tracer charge to unity. We also assume, without lack of generality, 
that E is oriented in the positive direction, i.e. > 0, and thus p> q. 



III. DEFINITIONS AND BASIC EQUATIONS. 

We start by writing the equations which describe dynamics of the tracer particle. Let 
XxiE^t) denote the position of the tracer at time t (by definition Xt{0) = 0) and P{X,t) 
be the probability that the tracer is at site X at time t. The mean displacement of the 



tracer particle, i.e. XT{E,t), is then defined as 



Xt{E) = ^XP(X,i) (4) 

X 

The time evolution of P{X, t) is governed by the equation 

P{X,t) = -P{X,t){prix+i{t) + qvx-i{t)) + 



+ rjxit) {p P{X - l,t) + qP{X+l,t)), (4) 
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where the dot denotes the time derivative. The first two terms on the r-h-s of Eq.(4) describe 
respectively the change in P{X,t) due to jumps of the tracer particle from the site X to 
sites X ±1, while the third and fourth terms account for jumps from the sites X ±1 to the 
site X. Multiplying both sides of Eq.(4) by X and summing over all lattice sites we arrive 
at the rate equation 

XT{E,t) = p/i - (5) 

where 

/, = J2Pi^^t)Vx+x{t) (6) 

X 

is the pairwise tracer-vacancy correlation function, which can be thought off as the prob- 
ability of finding at time t a vacancy at the distance A from the tracer. The correlation 
function is defined in the frame of reference moving with the tracer and jumps of the tracer 
change the value of fx. 

Consider now evolution of fx, which completely determines the mean displacement of 
the tracer, Eq.(5). The change in fx results from two different processes 

fx — Lbathifx) + Ltracifx), (7) 

where the operator L^ath accounts for the contribution coming from the motion of bath 
particles, whilst Lfrac describes the contribution of the tracer motion itself. Explicitly, for 

Lbathifx) we have (for |A| > 1) 

Uathifx) = - /a)(/a-1 + /a+i) - \fx{l - fx-l + 1 - /a+i), (8) 

where the first term describes the "birth" of a vacancy at the occupied site X + X due to the 
jumps of a bath particle from X + Xto vacant sites X + X±l; the second term describes a 
"death" of a vacancy at a vacant site X + X due to the jumps of a bath particle from sites 
X + X ± 1. One readily notices that the nonlinear terms in Eq.(8) cancel each other and 
Lbath is simply the second finite difference operator (for |A| > 1) 
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Uathifx) = ^{-2fx + fx-1 + fx+i) (9) 



The "diffusive" free particle like behavior in Eq.(9) is, of course, the consequence of the fact 
that all bath particles are identical, leading to cancellation of non - linear terms in Eq.(8). 

Consider now the contribution due to the motion of the tracer particle. In an explicit 
form we have (for |A| > 1) 



where the terms on the r-h-s of Eq.(lO) describe respectively the following events: (a) an 
occupied site at the distance A from the tracer becomes vacant if the tracer jumps into the 
previously vacant left-hand adjacent site and the site at distance A — 1 is vacant; (b) a vacant 
site at distance A from the tracer becomes occupied, if the tracer jumps into the previously 
vacant left-hand adjacent site and the site at distance A — 1 is occupied; (c) an occupied site 
at distance A becomes vacant when the tracer jumps into the previously vacant right-hand 
site and the site at distance A -|- 1 is vacant; eventually, (d) a vacant site at distance A 
becomes occupied when the tracer jumps into the previously vacant right-hand site and the 
site at distance A -I- 1 is occupied. 

Next, similar reasoning gives the behavior of f\ in the immediate neighborhood of tracer, 
i.e. at sites with A = ±1, which can be thought off as the boundary conditions for Eq.(7). 
Time evolution of f±i again is due to both the motion of the bath particles and of the tracer, 
i.e. can be represented in the form of Eq.(7), but here the operators Lj,ath and Ltrac are given 
as 




q /_i (1 - /,) fx-i - q /-I fx (1 - /a-i) + 



+ P /i (1 - fx) fx+i - P /i (1 - /a+i) fx, (10) 



Uath{f±l) = I f±2 (1 - /±l) - I f±l (1 - f±2), (11) 




P /i (1 - /2) + q f-i (1 - /i) (12) 



and 
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Ltrac{f-l) ^ Phil- /-l) - q /-I (1 - /-2) (13) 

Explicitly, we have for f±i the following equations 

/l = ^/2(l-/l) - ^/l(l-/2) - p/l(l-/2) + (14) 
/_! = i/_2(l-/-l) - ^/_l(l-/_2) + p/l(l-/-l) - ?/_i(l-/_2) (15) 

Another pair of boundary conditions for Eq.(lO) can be obtained by assuming that 
correlations between P{X,t) and r]x+x{t) vanish in the limit A ±oo, which seems quite 
plausible on physical grounds, and which may be checked numerically. We have from Eq.(6) 

/a^±oo = V±oo{t) E^(^'^) = Co, (16) 
X 

where we have used the normalization condition J2x — 1 aiid also the assumption 

that at infinitely large separations from the tracer the mean concentration of vacancies is 
unperturbed and equal to its equilibrium value Co- 

Eqs.(5), (7) to (13) constitute a closed system of equations for the derivation of Xx{E, t), 
within the framework outlined. Solutions to these equations will be discussed in the next 
section. 



IV. SOLUTION OF DYNAMIC EQUATIONS AND ANALYTICAL RESULTS. 

We first turn to the continuous-space description and rewrite our equations in the limit 
of continuous-A. Then, for A 7^ ±1 we have for the operator 

1 fx 

Lbathifx) ~ 2 'd)^' ^^'^'^ 



and, at points A = ±1, Eq.(ll) yields 

1 dfx 



Lbathif±l) 



, (18) 

A=±l 



2 dX 

Next, expanding /a_i ~ /a — dfx/ dX and /a+i ~ /a + dfx/ dX we obtain for the operator 
Ltrac, Eq.(lO), the following approximate expression 
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Ltracifx) -^f-^^ + Pf^^ = ^^(^'^) ^^^^ 

which holds for A 7^ ±1. 

The approximation as in Eq.(19), in which the coefficient before the gradient term is 
equal to the mean velocity of the tracer particle, is equivalent to the neglect of fluctu- 
ations in the tracer trajectory XT{E,t). In the following we show, however, that such 
a "deterministic" approximation is quite appropriate and leads to results which are in 
good agreement with simulation data. The reason this approximation works here is that 
in such a system the fluctuations of the tracer trajectories are essentially suppressed. 
This is due to the accumulation of the bath particles in front of the tracer which hin- 



ders fluctuations with XT{E,t) > XT{E,t). At the same time, the fluctuations with 



XxiE^t) < XT{E,t) are suppressed by the driving force. This leads to the mean-square 



deviation a^{t) = {XT{E,t)y — Xj^{E,t) growing only as (see the numerical results in 



Fig. 2, where the \Ja'^{t) is plotted versus in a striking contrast to the linear time de- 
pendence for conventional driven diffusive motion. Finally, within the same "deterministic" 
picture we obtain 



Ltrac{f±l) ~ TXT{E,t) (1 - /±i) (20) 

Combining Eqs.(17) to (20) we obtain the following continuous-space equation 

which is to be solved subject to the boundary conditions 

, 1 dh 



f. 



A=±l 



2 dX 



T XTiE,t) (1 - /±i) (22) 



A=±l 

and boundary conditions at A — > ±00, defined in Eq.(16). 



A. Solutions for A > 1. 

Let us consider first the behavior of /a for A > 1. We introduce a scaled variable 
u; — {X — 1)/Xt{E, t), (0 < a; < 00). In terms of u Eq.(21) takes the following form 
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"^'•^^ + a{u + l)^ = 0, (23) 



in which we have set 

d 



a = ^ {XT{E,t))\ (24) 
and which is to be solved subject to the boundary conditions 



fto=oo — Co, 



= a (1 - ^=o) (25) 

w=0 



We hasten to remark that such an approach presumes already that a is a time- 



independent constant and thus XT{E,t) grows as ^/at. However, such an assumption will 
be seen to be self-consistent if it turns out to be possible (and we set out to show that this 
is actually the case) to find the solution of Eq.(23) satisfying the boundary conditions in 
Eqs.(25). 

The appropriate solution to Eq.(23), which satisfies the boundary condition at infinity, 
reads 



= rdcu exp(-a(a; + a;V2)) + ^=o, (26) 

j+ Jo 



where 



poo I ^ I 

/+(«) = /+ = du; exp(-a(a; + o;V2)) = ^— exp(Q/2)[l - erf{^a/2)], (27) 

erf(x) being the error function. The value of /(^=o, which is the probability of having the 
r-h-s adjacent to the tracer site vacant, is to be defined from the boundary condition at 
(J — 0, Eq.(25), which gives 

Co — fu)=0 



1+ 



a{l - U=o) (28) 



We note here that a is as yet unspecified parameter. 



10 



B. Solutions for A < — 1. 



Consider next the behavior of fx, Eq.(21), past the tracer particle, i.e. in the domain 



A < —1. Now we choose the scaled variable 9 — —(A + 1)/Xt{E, t), {0 < < oo). In terms 
of 6 Eq.(21) takes the form 

^0. (29) 

with boundary conditions 

fe=oo = Co;^ /,=o) (30) 

aU \9=o 

The solution of Eq.(29), which satisfies the boundary at ^ = oo reads 
{Co - fe=o) 



fo = ^ ° , [ de exp(a {9 - 9^/2)) + /,=o, (31) 



where 



7_(a) = /_ = j^d9 exp(a {9 - 9'' /2)) = yjexp(Q;/2) [1 + erf{^f^)] (32) 

The constant of integration /g=o (which defines the probability of having the 1-h-s adjacent 
to the tracer site vacant) is to be determined from the boundary condition at ^ = 0, which 
gives 

^° ~ = - a (1 - fe^o) (33) 

Now we have a system of equations (28), (33) and (5), (24), which determine the behav- 



ior of XxiE^t), f^=Q and fe=o- In the following subsection we, construct a closed-form 
transcendental equation for the parameter a. 

C. Closed- form equation for a. 

We may formally rewrite Eq.(5) in terms of the parameters f^=o and f0=o as 



XT{E,t) = p^=o - qfe=^ (34) 
11 



Further, using Eqs.(34),(24) and (28) we have 

{Co — fu}=o) 



2I+XT{E,t) 



{P foj=0 - ?/e=o)(l - /a;=o) (35) 



Let us note now that the r.h.s. of Eq.(35) tends to zero as XT{E,t) grows. This means 
that (except for the case Co — 1, when Eq.(35) has only a trivial solution fu,=o — 1) the 
value pfoj=o should tend to qfe=o as t ^ oo. Expressing /^,=o as 

/,=o = ^ + 0, (36) 
P 

where is a deviation from the equilibrium value /a,=o, and assuming that at sufficiently large 
times this deviation is small enough to allow linearization of equation near the equilibrium, 
we have 

« - ^^^=° (37) 

2p7+(p - qfe=o) XT{E,t) 

Now, from Eq.(33) we express f0=o through the parameter a, which gives 

_ Co + al- , . 
1 + q; J_ 

and substituting Eq.(38) into the Eq.(37) we get 

^ go(p-g) + o^I-ipCo - q) 

2pl+{p - qCo + {p-q)aI_)XT{E,t) 



Finally, noticing that Xt {E,t) = p(j), and using the definition of the parameter a, Eq. (24) , 
we obtain the following closed transcendental equation for a 

Co(p-q) + a I- (p Co - q) 
" I+ip - qCo + {p-q)al^y ^ 

which, after some algebraic transformations, can be cast into a compact, symmetric form 

[a I+[a) + ) [a I-{a) + ) = — 41 

p-q p-q [p-qr 

where the a dependence of I± (Eqs. (27) and (32)) has been made exphcit. 

A simple analysis shows that Eq.(41) has positive bounded solutions a{p,Co) for any 

values of p, q — 1—p and Co (except the trivial case Co — 1 when bath particles are absent). 
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This means, in turn, that our approach to the solutions of coupled Eqs.(23),(25),(29) and 
(30) is self-consistent and XT{E,t) actually grows in proportion to y/t. In the trivial case 
when Co = 1 (when one expects, of course, the linear growth of XT{E,t)) Eq.(41) reduces 
to 

{a 1+ - 1) (a 7_ + 1) = (42) 

whose root is given by 

= (43) 

which means that a — oo and thus Xt{E, t) grows at a faster rate than y/t. 

Behavior of the prefactor a is presented in Fig.l, in which we plot the numerical solution 
of Eq.(41) for different values of the transition probabilities and different concentrations Cq. 
In the next subsection we present analytical estimates of the behavior of a in several limiting 
situations. 



D. Analytical estimates of limiting behaviors of the parameter a. 

We start by considering situations in which the parameter a can be expected to be small; 
namely, when p is close to q (in other words when p is only slightly above 1/2) or when the 
vacancy concentration is small, Co <^ 1. In the limit of small a the leading terms in Eqs.(27) 
and (32) behave as 



f -|- '-^ 



and Eq.(41) takes the form 



f¥a q-pCp fWa_ p-qCp _ p g (1 - Co)^ . . 
W 2 + p-g ^ W 2 + p-q > {p - qY ' ^^"^ 



and gives for small a, o; <C 1, 



" - -(%^)^ (46) 

TT i — Cq 
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Consequently, in this limit the explicit formula for the mean displacement will read 

^ C, (p- ,) /27 _ C„ tanh(^g/2) /2T ^^^^ 

1 — Co VTT 1 — Co VTT 

At this point it might be instructive to recall the analysis in which concerned the 
validity of the Einstein relation for one- dimensional hard-core lattice gases. Define the 
"diffusion coefficient" of the tracer particle in the symmetric case as -D = limt_^ooD{t), 



where D{t) = X^{E = 0,t)/{2t) and the "mobility" as /i = limE^oiU{E)/ E), where the 
stationary velocity is given by U{E) = limt^ooU{E,t), U{E,t) = XT{E,t)/t. In it was 
shown that for the infinitely large system the Einstein relation, i.e. the equation fi = j3D, 
holds trivially, since both D and /i are equal to zero. A somewhat stronger result has been 
obtained for the case when the one- dimensional lattice is a ring of length /. It was shown 
that here both D and fi are finite, vanish with the length of the ring as 1// and obey the 
Einstein relation = j3D{l). Our result in Eq.(47) reveals that for the infinite lattice 
Einstein relation is valid for the time-dependent "diffusion coefficient" and mobility, i.e. 
equation jj,{t) (= limE^o{U{E,t)/ E)) = [3D{t) holds exactly, provided that the time is 
sufficiently large such that both Eqs.(l) and (47) are valid. 

Consider next the situation when a is expected to be large, i.e. when C*o is close to unity. 
Expanding the functions /+ and /_ as 



1 1 



0, when a cxd, (4J 



and 



■exp(a/2) — > oo, when a oo, (49) 



a 



we arrive at the following equation 

(1 + - -) = — 7 — {V27ia exp(a/2) + ) (50) 

p — q a \P — q) P — q 

One may readily notice that when a is large, the r.h.s. of Eq.(50) tends to zero and thus we 
have 
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P(l-Co) ^ ' 

Eq.(51) allows formulation of the conditions when a is large more precisely: this happens 
when Co obeys the inequality 1 > Cq ^ q/p = exp(— 

In this limit the mean displacement of the tracer particle obeys 



XT{E,t) 



{p-q)t _ l{l-exp{-f3E))t ^^2) 



^p(l-Co) V I -Co 

Finally, let us consider the case of the directed walk in a one-dimensional lattice gas, 
when p = 1 and q = 0, which was previously examined in P^]. In this case Eq.(41) simplifies 
considerably and reads 

aoo/+ = Co (53) 
When Co is small, Cq <^ 1, we have from Eq.(53) 



too ~ 

TT 



(54) 



which is consistent with the result in Eq.(47), while in the limit Co ~ 1, Eq.(53) yields 

i — L-o 

i.e. exactly the behavior in Eq.(51) with p = 1 and q = 0. 



V. NUMERICAL SIMULATIONS. 

In order to check the analytical predictions of Sections IV and V Monte Carlo (MC) 
simulations have been performed. The simulation algorithm was defined as follows: The 
total number of particles in the system was held constant at 399 (398 bath particles and 
one tracer). We constructed a one- dimensional regular lattice of unit spacing and length 
2L + 1, sites of which were labelled by integers on the interval [—L,L]. For each Cp, the 
length was chosen as L = 7P. The number of particles and, thus, L was chosen to be large 
enough so that over the time scale of the simulations the mean vacancy density near L± was 
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unaffected by the dynamics of the tracer which always occupied the lattice site at i = 0. 
At the zero moment of MC time particles were placed randomly on the lattice with the 
prescribed mean concentration and the constraint that no two particles can simultaneously 
occupy the same site. The tracer particle was placed at the origin. The subsequent particle 
dynamics employed in simulations follows the definitions of Section 2 closely. At each MC 
time step we select, at random, one particle and let it choose a potential jump direction. If 
the selected particle is a bath particle, it chooses the direction of the jump - to the right or 
to the left, with the probability 1/2. If the selected particle is the tracer - it chooses the hop 
to the right with probability p and the hop to the left with probability q, q = 1 — p. For 
any particle the jump is instantly fulfilled if at this moment of time the adjacent site in the 
chosen direction is vacant. 

In simulations we followed the time evolution of the displacement of the tracer particle 
and plotted it versus the "physical time", which is the time needed for each particle to move 
once, on average. Fig. 2 presents typical behavior of the average displacement calculated for 
fixed q and Cq (the actual parameter for the presented results are q = and Cq = 1/3, the 
results for other values of Cq and q are similar) plotted versus y/i. The dotted line in this 
figure shows the evolution of the deviation from the mean average displacement Va^ which 
shows growth in proportion to t^^^. 

Numerical results for the time evolution of the mean displacement, performed for different 
values of parameters Cq and p, were used to obtain a numerical evaluation of the prefactor 
in Eq.(3) as a function of Cq and p. In Fig. 3 we plot, for the particular case of totally 
directed walk q — 0, the prefactor a versus the inverse concentration of vacancies. The solid 
line represents the results of our analytic calculations while the squares show the numerical 
results. In Fig. 4, for fixed concentration of the vacancies Cq = 1/3 we present the comparison 
between the analytical predictions of the prefactor a and the numerical MC results. 
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VI. CONCLUSIONS. 



To conclude, we have examined the behavior of the mean displacement of a driven 
tracer particle moving in a symmetric hard-core lattice gas. We have shown that the mean 
displacement grows in proportion to i.e. the hard-core interactions hinder the ballistic 
motion of the tracer and introduce effective frictional forces. The prefactors in the growth law 
are determined implicitly, in a form of the transcendental equation, for arbitrary magnitude 
of the driving force and arbitrary concentration of the lattice gas particles. In several 
asymptotic limits we find explicit formulae for these prefactors. Our analytical findings are 
in a very good agreement with the results of numerical simulations. 
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FIGURES 

FIG. 1. Theoretical values of the prefactor in the dependence of the average tracer displacement 

X{t) on t"*^/^ vs the probability of the "back step", q. 

FIG. 2. Typical dependence of squared mean tracer displacement and mean squared deviation 

on time. 

FIG. 3. Theoretical and experimental values of the prefactor in the dependence of the average 

tracer displacement, X{t), on i^/^ vs the inverse concentration of bath particles (q = 0). 

FIG. 4. Theoretical and experimental values of the prefactor in the dependence of the average 

tracer displacement, X(t), on t^/^ vs q {C = 1/3). 
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